IBM-1 calculations towards the neutron-rich nucleus 106 Zr 
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The neutron-rich N = 66 isotonic and A = 106 isobaric chains, covering regions with varying types 
of collectivity, are interpreted in the framework of the interacting boson model. Level energies and 
electric quadrupole transition probabilities are compared with available experimental information. 
The calculations for the known nuclei in the two chains are extrapolated towards the neutron-rich 
nucleus 106 Zr. 

PACS numbers: 



I. INTRODUCTION 

In the last decade the neutron-rich nuclei in the 40 < 
Z < 50 region have attracted both theoretical and 
experimental attention. They were extensively stud- 
ied via spontaneous or induced fission reactions. Nu- 
clei from this region of Segre chart exhibit vibrational, 
transitional, and rotational types of collectivity. The 
neutron-rich palladium and ruthenium nuclei, for exam- 
ple, show a typical transitional behavior while molyb- 
denum isotopes exhibit a vibrational to near-rotational 
evolution Such changes in the degree of collectiv- 

ity are even stronger in the zirconium isotopic chain [l|, Q 
where the structure evolves from near- vibrational in 98 Zr 
to rotational-like in 104 Zr, the latter being the most de- 
formed nucleus of the neutron-rich Z = 40 isotopes. As 
the next zirconium isotope 106 Zr is an N — 66 mid-shell 
nucleus, an even larger deformation can be expected. 

A detailed analysis of some spectroscopic observables, 
such as the ratio R4/2 of excitation energies of the first 2 + 
and 4+ levels fli] or the amplitude of even-odd stagger- 
ing in the 7 band , shows that mid-shell effects in pal- 
ladium isotopes occur in 114 ' 116 Pd, i.e. two and four neu- 
trons beyond N = 66. In the ruthenium isotopes these 
effects arise precisely at N = 66 while in the molybde- 
num chain they occur in 106 Mo, i.e. two neutrons before 
mid-shell. Furthermore, it was shown that the degree of 
collectivity decreases in the heavier N = 66, 68 molyb- 
denum isotopes 108 ' 110 Mo. In spite of the fact that the 
experimental information on the heavy zirconium nuclei 
is poor, it is clear that the degree of collectivity increases 
towards 104 Zr [l[ which has 64 neutrons and is the heav- 
iest zirconium isotope known to date. The shift of the 
mid-shell in the various isotopic chains was explained in 
Ref. m as resulting from the filling of the neutron h\i/2 
and proton g 9 / 2 intruder orbitals. 

Self-consistent calculations show that in the region 
of neutron-rich nuclei approaching the neutron-drip 
line, the single-particle shell structure may significantly 
change due to the diffuseness of the neutron density [3| . 
Also, as widely discussed in the literature, due to a pre- 
dicted weakening of the spin-orbit force, new magic num- 
bers may be expected in neutron-rich nuclei. For in- 



stance, in a study of N = 82 isotones, evidence for shell 
quenching was presented in Ref. [J]. On the other hand, 
Jungclaus et al. [f| argued that from the isomeric decay 
study in 130 Cd there is no direct evidence for the shell 
quenching in the region of the heavy cadmium nuclei. 

The lower-Z N = 82 nuclei, however, have not yet 
been studied experimentally and the suppression of the 
shell effects in this region is still an open question. Prop- 
erties of these nuclei were predicted within a relativistic 
Hartree-Fock Bogoliubov (HFB) approach Q. The calcu- 
lations, performed with different parametrizations, show 
that the N = 82 shell gap persists in the heavy palla- 
dium and ruthenium nuclei but that a weakening of this 
gap is expected in the zirconium isotopes. For example, 
the two-neutron separation energies from 118 Zr to 126 Zr, 
obtained via HFB calculations with the SkP interaction, 
show no discontinuity at N = 82 but a rather smooth be- 
havior. As discussed above, nuclear collective behavior is 
expected to enhance towards the mid-shell as a result of 
the increase in valence particle number. If, however, the 
HFB+SkP scenario for the zirconium isotopes turns out 
to be valid and N — 82 vanishes or is weakened in the 
neutron-rich region, then the determination of valence- 
particle number is fraught with ambiguity, leading to a 
different behavior of nuclear collectivity in this isotopic 
chain. The aim of the present paper, therefore, is to pre- 
dict the spectroscopic properties concerning the N = 66 
mid-shell zirconium nucleus 106 Zr, based on the assump- 
tion that N — 82 remains a magic number in the heavy 
zirconium isotopes. Comparison with results of future 
experiments on this nucleus will then reveal whether this 
hypothesis is borne out or not. 



II. MODEL 

Often the only available information about neutron- 
rich nuclei comes from prompt 7-ray spectroscopy. The 
observed levels are grouped into A J = 2 sequences, corre- 
sponding to the ground-state band, and/or into AJ = 1 
sequences, based on a l\ level which usually is inter- 
preted as a quasi-7 band Q. Such structures naturally 
appear in the framework of the interacting boson model 
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(IBM) 8] which has been shown to be successful in the 
description of nuclear collective properties. 

The IBM in its first version, known as IBM-1, is based 
on the assumption that nuclear collectivity can be ex- 
pressed in terms of s and d bosons [8| . The model Hamil- 
tonian is constructed from a set of 36 operators, bilin- 
ear in the boson creation and annihilation operators and 
generating the U(6) Lie algebra. Dynamical symmetries 
occur if the Hamiltonian can be written as a combina- 
tion of invariant (or Casimir) operators of specific sub- 
algebras of U(6) and three such cases occur, namely 
the spherical vibrational limit U(5), the deformed limit 
SU(3), and 7-soft limit SO (6). These dynamical symme- 
tries generate energy spectra with states that are labeled 
by the irreducible representations of the algebras in the 
respective chains that reduce the dynamical algebra U(6) 
into the symmetry algebra SO (3) of rotations in three- 
dimensional space. The different limits thus correspond 
to nuclei with distinct collective properties. For example, 
in the SO (6) limit the ground-state and 7 bands lie in the 
same representation while in the SU(3) limit they are in 
different ones, leading to forbidden or weaker inter-band 
transitions in the latter case. 

In the A ss 110 nuclei of interest here strong transitions 
between the two bands are observed indicating that none 
of these nuclei can be interpreted in the exact SU(3) limit, 
but rather a transitional behavior should be expected. 
This can be achieved in the IBM-1 by the use of the full 
Hamiltonian which reads Q 

H = edfid + kQ x -Q x + k'L ■ L + C3T3 • T 3 + C4T4 • T 4 , (1) 

where fid = eft - J and L M = \/l0[^ x dyp are the <i-boson 
number and the angular momentum operators, respec- 
tively. Furthermore, the quadrupole operator is defined 
as Q* = [eft x 5]^ +x[(ft x d\ii 2 \ while the last two terms 
in the Hamiltonian involve the operators = [(ft x d]^ 3 ' 
and Ti il± = [eft x d]^ . In the present work there is no 
need for the rotational term L ■ L and we take k' = 
throughout. The total number of bosons N = n s + rid 
is taken as half the number of valence particles or holes, 
counted from the nearest closed-shell configuration, fol- 
lowing the prescription of Ref. [l(|. As such, IBM-1 cal- 
culations are indirectly related to the underlying shell 
structure. Finally, in the consistent-Q formalism [ll| the 
operator for electric quadrupole transitions is introduced 
as T^(E2) = e^Q x , where eb is the boson effective charge 
and Q x is the quadrupole operator with the same param- 
eter x as in the Hamiltonian fl}. 



III. RESULTS AND DISCUSSION 

Since the addition or subtraction of a few nuclcons 
may change significantly a nuclear spectrum, two nuclei 
with different numbers of valence neutrons and protons 
but with the same total number of valence nucleons — 



TABLE I: Parameters and rms deviation (in keV) for the 
N = 66 isotones and the A — 106 isobars. 
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-0.30 


58 


A = 


106 


1053 


— 25.6 


-70.6 


-118.3 


-0.30 


143 
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and hence the same number of bosons — may display dif- 
ferent types of collectivity. In adjusting the parameters 
of the Hamiltonian ([1]) to spectra observed in a given 
region of the nuclear chart, it is therefore not sufficient 
to assume them to be constant for all nuclei since that 
would lead to identical spectra for nuclei with the same 
N. Instead, some dependence of the parameters on the 
separate neutron and proton boson numbers N v and N„ 
should be imposed. To avoid the complexity of deter- 
mining the correct functional dependence on N v and N n 
for each of the Hamiltonian parameters, we have followed 
the simpler procedure of separately considering chains of 
isotopes, isotones, or isobars that end with the nucleus 
of interest, 106 Zr, and fitting each chain with an indepen- 
dent set of parameters. This method avoids the problem 
of having identical spectra for constant boson number N 
(such nuclei belong to different chains) and, in addition, 
since the extrapolation to 106 Zr obtained from the differ- 
ent chains may vary, it gives an idea on the possible error 
of the prediction. 

In the present application to 106 Zr, the N = 66 isotonic 
chain includes 7-unstable 112 Pd, more triaxial 110 Ru, 
and near-rotational 108 Mo. The second set of nuclei be- 
longs to the A — 106 isobaric chain [l6| extending from 
vibrational-like 106 Pd to near-rotational 106 Mo. Unfor- 
tunately, the zirconium isotopic chain does not allow the 
determination of a unique set of model parameters since 
only yrast data are available for the neutron-rich isotopes 
100,102 ' 104 Zr. Consequently, only the isotonic and isobaric 
chains are considered in the fit. Nuclei in a given chain 
are distinguished by the number of bosons (i.e., parti- 
cle or hole pairs), counted from the nearest closed-shell 
configurations which are Z — 28 or Z — 50 for the pro- 
tons, and N = 50 or N — 82 for the neutrons. Figure [1] 
summarizes known nuclei in the neighbourhood of 106 Zr. 
Nuclei with known excited levels are distinguished (full 
lines) from those of which only ground-state properties 
are known (dashed lines). Furthermore, the nuclei used 
in the fit are high- lighted with thick lines. 

The model parameters ed, k, C3, and C4, listed in Ta- 
ble U are obtained after a least-squares fit to the exper- 
imental level energies. Since its effect on the excitation 
energies is weak, the parameter \ is varied independently 
to give a fair description of E2 transition rates in 106 Pd 
within the consistent-Q formalism, and other nuclei are 
calculated with the same value. Although within a given 
chain parameters may vary with boson number, any al- 
lowed variation does not markedly improve the quality of 
the fit. Consequently, to keep the number of parameters 
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FIG. 1: Even-even nuclei in the region close to Zr. All nuclei in boxes with full lines have known excited levels while those 
in boxes with broken lines are only known in their ground state. The nuclei used in the fit are in boxes with thick lines. 



to a minimum, they are taken constant for a given chain 
and any structural evolution within that chain is due to 
the changing boson number N. 

The quality of the fit is summarized by the root-mean- 
square (rms) deviation between the experimental and cal- 
culated level energies. The rms values obtained in the two 
independent fits are given in Table HI The experimen- 
tal mm in mm 

and theoretical energies are listed 
in Tables III] and HIT] for the N = 66 and A = 106 chains, 
respectively. The nucleus 106 Pd is known in greater de- 
tail and is shown separately in Table IIVI In both the 
isotopic and isobaric chains, the observed energy of the 
2f state decreases when approaching the nucleus 106 Zr. 
The same trend is observed for the calculated energies 
as a function of boson number. The experimental en- 
ergy of the 7-band head also decreases with increasing 
boson number and a similar evolution is observed for the 
theoretical energies. According to the systematics of 
the less neutron-rich nuclei, the 2^ level in 106 Zr should 
thus lie around 140 keV while the 2^~ level should be in 
the region 600-650 keV. Furthermore, the two different 
fits give approximately the same 7-band energies. Even 
though the two fitted parameter sets are close (see Ta- 
ble |TJ), the N = 66 set gives rise to a slightly faster in- 
crease of the ground-state band energies than is the case 
for the A = 106 set. 

The displacement of odd-spin with respect to even- 
spin 7-band levels can be quantified with the following 
function: 

S(J) = £(J) _(i±iW__Mi±i), (2) 

where E(J) is the energy of the level with angular mo- 
mentum J. The experimental and theoretical staggering 
in the N = 66 isotones is plotted in Fig. [2j to gauge any 
systematic behavior, plots are drawn to the same scale. 
In general, the staggering amplitude is high for spherical 
nuclei and decreases towards the mid-shell region. The 



N = 66 isotones do not constitute an exception to the 
systematic trends established in Ref. [|J • The 7-unstable 
nucleus 112 Pd (10 bosons) shows the highest amplitude 
in S(J) among the N = 66 isotones. The 7-band stag- 
gering is less pronounced in 110 Ru (11 bosons) while in 
108 Mo (12 bosons) it is highly suppressed. From these 
systematics trends a 7 band with a very low staggering 
amplitude is expected in 106 Zr, namely S(J) ~ 10 keV, 
which would make 106 Zr the nucleus with the lowest stag- 
gering amplitude observed in the 40 < Z < 50 region. 

Also electric quadrupole transition probabilities have 
been considered for the neutron-rich nuclei belonging to 
the N = 66 and A = 106 chains. The calculated val- 
ues are listed in Tables [V] and IVII along with the few 
known experimental values (l6j . The boson effective 
charge eb = 0.11 eb is obtained after a least-squares fit 
to the B(E2) values measured in 106 Pd and this value is 
used for all other nuclei. The results show an increase 
of the £(E2;2+ -> 0+) value from the near- vibrational 
106 Pd towards the near- rotational 108 Mo. This fl(E2) 
value is predicted highest for 106 Zr and the calcula- 
tions for the two different chains give a consistent value 
of £(£2; 2+ -> 0+) fa 0.50-0.55 e 2 b 2 . Furthermore, 
we note from Table IVII the behavior predicted for the 
B{E2; 2+ -> 0^) value which is small in 106 Pd and 106 Zr 
but becomes large for the isotopes in between, a charac- 
teristic feature of transitional nuclei. 



IV. CONCLUSIONS 

The aim of this paper was to illustrate the usefulness of 
the interacting boson model in its simplest version, the 
IBM-1, for predicting properties of exotic nuclei. The 
usual difficulty that arises with such attempts is that 
the variations of the model's parameters with the va- 
lence neutron and proton numbers are not known, pre- 
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FIG. 2: Experimental (full symbols connected with full lines) and theoretical (open symbols connected with dashed lines) 
7-band staggering S(J) in the N = 66 isotones. In 106 Zr both near-identical predictions are shown based on the extrapolation 
of the N = 66 isotones (open squares) and of the A = 106 isobars (open diamonds) . 



TABLE II: Experimental and theoretical excitation energies 
(in keV) of levels in neutron-rich N = 66 isotones. 



nucleus 




T-iex 

Ksb 


gith 




Z? ex 


TTlth 

7 




o+ 


o 


o 


2+ 


736 


775 


(N = 10) 


2+ 


348 


354 


3+ 


1096 


1138 




4+ 


883 


849 


4+ 


1362 


1304 




6+ 


1550 


1481 


5+ 


1759 


1692 




8+ 


2318 


2240 


6+ 


2002 


1934 




10+ 


3049 


3118 


7+ 
8+ 
9+ 
10+ 


2483 
2638 
3085 
3327 


2357 
2669 
3129 
3507 


iil) Ru 


0+ 








2+ 


613 


689 


(iV = ll) 


2+ 


241 


253 


3+ 


860 


960 




4+ 


663 


671 


4+ 


1084 


1139 




6+ 


1239 


1239 


5+ 


1375 


1449 




8+ 


1945 


1945 


6+ 


1684 


1699 




10+ 


2759 


2777 


7+ 
8+ 
9+ 
10+ 


2021 
2397 
2777 
3255 


2060 
2377 
2789 
3168 


u)8 Mo 


0+ 








2+ 


586 


647 


(N = 12) 


2+ 


193 


178 


3+ 


783 


830 




4+ 


564 


533 


4+ 


978 


1010 




6+ 


1090 


1047 


5+ 


1232 


1257 




8+ 


1753 


1705 


6+ 


1508 


1505 




10+ 


2529 


2495 


7+ 
8+ 
9+ 
10+ 


1817 
2170 
2524 
2950 


1819 
2132 
2506 
2880 


lU6 Zr 


0+ 








2+ 




618 


(N = 13) 


2+ 
4+ 




141 

455 


3+ 
4+ 




754 
922 




6+ 




926 


5+ 




1133 




8+ 




1543 


6+ 




1371 




10+ 




2297 


7+ 
8+ 
9+ 
10+ 




1653 
1956 
2303 
2668 



TABLE III: Experimental and theoretical excitation energies 
(in keV) of levels in neutron-rich A = 106 isobars. 



nucleus 


V 


-C'Rsb 


TTlth 
-C'Ksb 




77. ex 


7-ith 


il)b Ru 


0+ 








2+ 


792 


772 


(AT = 9) 


2+ 


270 


362 


3+ 


1092 


1156 




4+ 


715 


865 


4+ 


1307 


1308 




6+ 


1296 


1498 


5+ 


1641 


1720 




8+ 


1973 


2252 


6+ 


1908 


1943 




10+ 


2705 


3115 


7+ 
8+ 


2284 
2960 


2387 
2676 


i,lb Mo 


0+ 








2+ 


711 


651 


(AT =11) 


2+ 


172 


205 


3+ 


885 


881 




4+ 


522 


582 


4+ 


1068 


1057 




6+ 


1033 


1110 


5+ 


1307 


1338 




8+ 


1688 


1773 


6+ 


1563 


1578 




10+ 


2472 


2559 


7+ 
8+ 
9+ 
10+ 


1868 
2194 
2559 
2951 


1916 
2217 
2608 
2965 


iub Zr 


0+ 








2+ 




630 


(N = 13) 


2+ 
4+ 




134 
433 


3+ 
4+ 




762 
924 




6+ 




883 


5+ 




1129 




8+ 




1473 


6+ 




1357 




10+ 




2192 


7+ 
8+ 
9+ 
10+ 




1630 
1919 
2254 
2601 



eluding a reliable extrapolation to unknown regions of 
the nuclear chart. Our method proposes to circumvent 
this problem by studying the structural evolution in three 
different chains of nuclei, namely isotopic, isotonic, and 
isobaric ones, which cross at the nucleus of interest at 
the outskirts of the region of stable nuclei. Predictions 
are obtained by extrapolating the different chains to the 
exotic nucleus in question, and, in addition, a compar- 
ison of these extrapolations gives an idea of the errors 
involved. 
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TABLE IV: Experimental and theoretical excitation energies 
(in keV) of levels in 106 Pd. 



nucleus 


r 


E 1 


E2 


E3 


£4 


E 5 


iu« pd 


o£ 





1134 


1706 


2001 


2278 


(N=7) 


o t + h 





1017 


1384 


2080 


2297 




2+ 


512 


1128 


1562 


1909 


2242 




2+ 


519 


982 


1575 


1859 


2004 




3+ 

°ex 


1558 














1483 












4+ 


1229 


1932 


2077 


2283 






4+ 


1153 


1617 


2010 


2213 






5i 


2366 












5+ 

°th 


2141 












6+ 
6 t + h 


2077 
1896 











TABLE V: Theoretical B(E2) values (in units e 2 b 2 ) for tran- 
sitions in neutron-rich N = 66 isotones. 



transition 


112p d 


11() Ru 


108 Mo 


ioe Zr 


2+^0+ 


0.27 


0.35 


0.45 


0.53 


4+^2+ 


0.43 


0.54 


0.65 


0.75 


0+^2+ 


0.22 


0.17 


0.057 


0.003 


2+^0+ 


0.015 


0.024 


0.025 


0.020 



This method was applied to the neutron-rich members 
of the N = 66 isotonic and A = 106 isobaric chains of 
which level energies and electric quadrupole transition 
probabilities were fitted with IBM-1. The two chains in- 



tersect at 106 Zr which allowed the prediction of this nu- 
cleus' excitation energies and electric quadrupole transi- 
tion properties. 

As a final remark we emphasize that the IBM-1 is a 
valence-nucleon model and that the extrapolations as de- 
scribed here crucially depend on the definition of neigh- 
boring closed-shell configurations. Thus, the results on 
106 Zr indirectly involve the assumption of the persistence 
of magic numbers in this region of neutron-rich nuclei. 
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TABLE VI: Experimental and theoretical B(E2) values (in 
units e 2 b 2 ) for transitions in neutron-rich A — 106 isobars. 



transition 


106 pj 




106 Ru 


106 Mo 


106 Zr 




Expt 


Theo 


Theo 


Theo 


Theo 


2+-0+ 


0.134 (4) 


0.14 


0.23 


0.37 


0.54 


4+^2+ 


0.211 (17) 


0.22 


0.36 


0.55 


0.77 


0^-2+ 


0.105 (20) 


0.13 


0.16 
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0.0044 


2+^0+ 


0.0035 (3) 


0.0019 


0.0097 


0.023 


0.021 
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